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Abstract 

The electrostatic force exerted on a targeting sequence by the electrical potential across the 
inner mitochondrial membrane is calculated and found to vary from 1.4 pN to 2.2 pN (per unit 
elementary charge) as the radius of the inner membrane pore (assumed aqueous) is varied from 
12 to 6.5 A, its measured range. Since the pore is not very much wider than the distance between 
water molecules, the full shielding effect of water may not be present; the extreme case of a non- 
aqueous pore gives a force of 3.1 pN per unit charge, which represents an upper limit. When 
applied to mitochondrial import experiments on the protein barnase, these results imply that a 
force of 11 ± 4 pN is sufficient to catalyze the unfolding of barnase during import. Comparison 
of these results with unfolding forces measured using atomic force microscopy suggests that the 
two are not inconsistent. 


1 Introduction 

Most mitochondrial proteins are encoded in the cell’s nuclear DNA, manufactured within the cytosol 
as precursors, and translocated into mitochondria across the organelle’s inner and outer membranes 
(henceforth IM and OM) through points where the two membranes come into contact. For the 
basic facts consult Alberts et al. (1994) and Pfanner and Neupert (1990). After translocation, the 
precursors are sent to the appropriate mitochondrial subcompartment where they are assembled into 
protein complexes. Most precursors that are targeted to the lumen of the mitochondria, called the 
matrix, are synthesized with a targeting sequence (TS), also called a presequence, attached at their 
amino terminus. This TS marks the precursor for translocation. We are concerned with precursors 
that are folded prior to import and where the TS protrudes from the precursor. Targeting sequences 
of this kind always have an abundant number of positively charged residues with few negative ones. 
As previously suggested (e.g , Martin, Mahlke, and Pfanner, 1991) the positive charges allow the 
inner membrane’s electric potential to exert a force that is directed into the mitochondrion. 

The translocation of protein precursors into mitochondria involve a number of actors (Pfanner 
and Truscott, 2002) besides the membrane potential; see Fig. 1. The TS first interacts with protein 
receptors (Tom20 and Tom22) on the surface of the outer membrane. These receptors may promote 
insertion of the TS into the OM pore, which itself consists of the protein Tom40. The pore of the 
inner membrane likewise consists of transmembrane proteins (Timl7 and Tim23). A portion of the 
Tim23 protein that lies exposed on the outer face of the IM appears to facilitate insertion of the TS 
into the IM pore; the membrane potential activates the insertion (Bauer et al., 1996). The passage 
of the TS through the IM pore may be driven by thermal motion, the electric field of the membrane 
potential, interaction with the Tim proteins, or a combination. 
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Figure 1: Schematic of protein import. 


Once the TS has been threaded into both OM and IM pores, the bulk of the protein lying 
on the outer mitochondrial surface must then unfold. Huang et al. (1999) concluded that the 
unfolding is initiated at the targeting sequence and that precursor proteins are unraveled sequentially 
from their N- ter mini. The unraveling occurs when the targeting sequence engages the unfolding 
machinery associated with the inner mitochondrial membrane while the structured domain remains 
at the entrance to the import channel. The simplest mechanism by which the import machinery 
could unravel a protein at a distance would be by pulling at the targeting sequence. Atomic force 
microscopy (AFM) experiments show that the N-terminus of a protein needs to be pulled only a 
short distance before the protein denatures. This distance is an empirically defined width of the 
potential well for unfolding and its values range between 3 and 17 A for different domains (Rief et 
al., 1997; Rief et al., 1998; Best et al., 2001). 

What pulls the targeting sequence through the required distance? If the TS is long enough 
to span both membranes and reach sufficiently far into the mitochondrial matrix, then Tim44 in 
association with mtHsp70 is able to unfold the protein by an ATP driven action (e.g., Matouschek, 
Pfanner, and Voos, 2000). Many targeting sequences, however, are not long enough to span both 
membranes; for instance the total thickness of yeast mitochondrial membranes is at least « 140 A. 
This corresponds to 40 amino acids in the fully extended conformation while the average length of 
yeast presequences is smaller, viz, 31 amino acids (Huang, Ratliff, and Matouschek, 2002). When 
targeting sequences are not long enough to interact with mtHsp70, the rate of import of precursor 
proteins depends upon the strength of the electrical potential and the number of positively charged 
amino acids (Huang, Ratliff, and Matouschek, 2002). The simplest implication of this result is that 
for short targeting sequences, the force exerted by the inner membrane potential upon the charged 
residues of the targeting sequence unfolds the passenger protein. In this work we investigate this 
hypothesis by calculating the electrostatic force exerted by the potential and make a preliminary 
attempt to determine whether it is sufficient to unravel a protein. 
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Figure 2: Sketch for computational model. Subscripts: o, outer membrane; i, inner membrane; c, 
charge layer. 


2 Model 

2.1 Computational Model 

A distribution of charge density p(x) (per unit volume) in a medium with dielectric constant e(x) 
produces an electric field E — — V0, where 0 is obtained from 

V . (cV0) = -4 t rp. (1) 

Fig. 2 shows the geometry assumed for the distributions of e(x) and p(x). Shading denotes mem- 
braneous regions where e = 1. Lack of shading denotes aqueous regions (e = 80). These include the 
cytosol , mitochondrial matrix, and the intermembrane space where the layer of + charge is located. 
The dielectric constant incorporates the shielding of charges by water dipoles. There is also screening 
by ions; at pH near 7, however, the Debye shielding distance (e.g., Probstein, 1994) is about 1 /im, 
rendering this effect negligible at the length scales of present interest. 

The OM and IM pores are also assumed to be aqueous based on their observed hydrophilic 
character (see Hill et al., 1998; Truscott et al., 2001). The radius r Q of the OM pore was taken to 
be 12 A based on reported measurements, namely, between 10 and 13 A according to Schwartz and 
Matouschek (1999), 11 A according to Hill et al. (1998), and 10 A according to Kunkele et al. (1998). 
Less is known about the radius r* of the IM pore. Schwartz and Matouschek (1999) concluded that 
n is at most 10 A. Here we will consider values in the wider range 6.5 < r* < 12 A, suggested by 
Truscott et al. (2001). Since the spacing between water molecules is about 3 A, only a few water 
molecules will be able to occupy the pores. To qualitatively allow for such an effect, values of 
down to 0 A, representative of a non-dielectric pore, will also be considered. 

For the value of the mitochondrial membrane potential, A#, we used 150 mV corresponding to 
the protein import experiments of Huang, Ratliff, and Matouschek (2002). The charge density a 
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8 = 1 

Figure 3: Sketch for analytical model. 


(per unit area) is then inferred to be 


A# 

“ 4tt(2g0 ’ 


( 2 ) 


where 2 d = hi + h c is the distance between the charged layers. The volume charge density is then 

.p — a jh c 

We took hi = h 0 = 65 A, a value consistent with electron micrograph pictures. The radius r c of 
the holes in the charged layers was taken to be 10 A. The thickness of charge layers h c = 10 A, a 
reasonable guess; the electric, field away from the edges of the charge layers will be insensitive to h c . 

Equation 1 was solved numerically using a B-spline Galerkin scheme (Shariff and Moser , 1998) 
in cylindrical polar coordinates (x,r), where x is the axial coordinate (measured from the entrance 
of the IM pore and positive into the mitochondrion) and r is the radius. The discretization cells 
were designed to be small at interfaces where jumps in dielectric constant and charge density occur, 
and to become larger as the computational boundary is approached. In most runs the smallest 
computational cell size was lA x lA and the computational domain was x E [—200, 200] A, r E 
[0, 110] A. As a check on accuracy, a computation with half the cell sizes in each direction and twice 
the radial domain size was also run. The boundary condition d<f>/ dn — 0 was applied at the boundary 
of the computational domain which is large enough for the boundary condition to be accurate. Here 
n is the coordinate normal to the boundary. At the symmetry axis we also required d<p/dr = 0 which 
is precisely the condition required for an axisymmetric function to have continuous radial derivatives 
at the axis. Since the Galerkin method is based upon integrals, discontinuous distributions of e(x) 
and p(x), which occur in the present model, can be treated. At an interface across which e suffers a 
jump, E n , the component of the electric field normal to the interface also jumps. Since the computed 
solution is a projection of the exact solution upon the space of B-splines, this jump leads to some 
Gibbs oscillation in E n . Such oscillation may be witnessed in Fig. 5 and was generally found to be 
weak. 


2.2 Analytical Model 

Since V x E = 0, the tangential component of E is always continuous across charge layers and across 
discontinuities in e; since V (cE) = 0 outside of charge layers, the normal component of E suffers a 
jump across discontinuities in e (see e.g., Jackson, 1962, Ch. 1). In particular, when e increases by a 
factor of 80 across an interface, the component of the electric field normal to the interface diminishes 
by the same factor. 

For the case of a thin pore (r* hi and r 0 <& h Q , which is typical), we expect that the electric 
field in the membrane space, will be primarily in the axial (x) direction, and will therefore persist 
in the aqueous pore without being substantially diminished by the presence of water. Hence, as an 
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approximation we take t — I. everywhere. In addition we take the charge layers to be infinitesimally 
thin plates (Fig. 3). _ 

In this sub-section, let the positive axial direction (£) point out of the mitochondrion. Let E(£) 
denote the axial (£) component of the electric field due to a single plate with a hole (at £ = 0) of 
radius r c and surface charge density a. We are evaluating the electric field along the axis of the hole 
(r = 0) where only the axial component is non-zero. Let E'(£) denote the axial electric field of the 
charged disk which closes the hole. Then, 

£(0 + £'(0 = 2tt<7, (3) 


the right-hand side being the electric field of an infinite plate. Integrating the inverse square law 
over the disk gives 

where cos# = £/\/£ 2 + r ^* Carrying out the integration in (4) and using (3) gives 

27r 


£(**) = /“TT72- 

V r c + r 

Finally, superposing two plates spaced a distance 2d apart we have: 

E(x) = E(x — d; a) + E(x + d; — cr), 


(5) 


( 6 ) 


where x is measured from the midpoint of the two charged layers (see Fig. 3). The peak value of 
the field occurs at x = 0 anc. is: 

Er\ 


inad 


" peak vVc 2 + d- ' 


^npak — 47T<T - 


(7) 

(8) 


In the limit r c « d we get: 

Vak - ’ 2d ’ 

the electric field in a parallel plate capacitor. The corresponding axial force F eiPea k = e£ pea k P er 
unit (+) elementary charge : e) is 


-^e,peak(P-^) 


- 1.6 


A$(mV) 

2d(A) ’ 


(9) 


in the ^ direction. 


Equation 9 is a convenient formula for calculating an upper limit on the force. 


3 Results and Discussion 

Fig. 4 is a result of the computational model. It shows the force field exerted on a particle having 
a charge equal to that of the electron but with a positive sign. Colors depict the magnitude of the 
force and the arrows provide its direction. The force within the IM pore is remarkably uniform, 
both radially and axially. Whatever leakage there is of the field at the entrance of the IM pore has 
a direction that is favorable o centering and insertion of the targeting sequence into pore. The field 
direction at the exit of the IM pore is favorable for diffusion and exit out of the pore. 

The radial uniformity of the force field within the pore (r < r* = 6.5 A here) is illustrated in 
Fig. 5. At the pore boundary (r = 6.5 A) the force field suffers a jump in derivative and then at 
large distances from the pore relaxes slowly to the uniform field of a parallel plate capacitor. 

Fig. 6 plots the force along the axis of the pores for various values of the radius of the IM pore. 
As the radius, r t , of the IM pore decreases, the electric field increases as a result of less shielding by 
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Figure 4: Electrostatic force per elementary (+) charge. The radius, r* of the IM pore is 6.5 A here. 
Colors depict the magnitude of the force while arrows show its direction. 
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Figure 5: Radial profile of the force per elementary charge for the case of r* = 6.5 A. The profile is 
taken at the mid-section of t he pore. The small oscillation is Gibbs phenomenon, an artifact of the 
numerical solution. 


water. The result (solid line) of the analytical model (6) provides an upper bound and becomes a 
better approximation to the results of the computational model as the pore radius diminishes. The 
force profile given by the computational model is quite uniformly distributed along the pore and 
indicates a small leakage of the electric field into the OM pore. If we let the radius of the IM pore 
take on values ranging from 0 A (to allow for the possibility of a completely non-dielectric pore) up 
to 12 A, then we conclude that 1.4 < F ey pea k <3.1 pN. 

Next, consider import experiments for the protein barnase (Huang, Ratliff, and Matouschek, 
2002) and focus on those corses in which the unfolding mode ranges from being spontaneous to 
membrane potential driven (see Table 1), i.e., exclude cases for long presequences which are unfolded 
by mtHsp70. The column in Table 1 labeled “Unfolding mode” gives a qualitative indication of 
the relative importance of spontaneous (denoted S) and membrane potential (denoted M) driven 
unfolding in each experiment . The net number of positive charges is determined as described in the 


Targeting 

Sequence 

Net no. of + charges 
in IM pore 

Total force 
(pN) 

Import rate 

(A s- 1 ) 

Unfolding mode 

(35) 

3 

4.2 9.3 

0.35 ± .02 

ssss 

(35; A16K) 

4 

5.6-12.4 

0.53 ± .04 

SSMM 

(35+5) 

4 

5.6-12.4 

0.87 ± .09 

SSMM 

(65) 

0-3 

0.0-9. 3 

1.3 ± .2 

SMMM 

(35; E15L) 

5 

7.0-15.5 

2.4 ± .1 

MMMM 


Table 1: Predicted forces fo mitochondrial import experiments (Huang, Ratliff, and Matouschek, 
2002) 
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Figure 6: Force per elementary (+) charge. The origin of the abscissa lies at the entrance of the 

IM pore which is 65 A long. Positive force is directed into the mitochondrion. : using the 

analytical result (6); computational results: , r* = 3 A; , r* = 6.5 A; , rq = 12 A. 


appendix and after multiplication by the range of F e in the previous paragraph, we obtain the values 
in the “Total Force” column. One observes a trend, more or less, of increasing import rate with 
increasing predicted force; the case of targeting sequence (65) is anomalous. Import with the (35; 
E15L) presequence is, by all indications of the experiment, catalyzed by the membrane potential; 
the force that accomplishes this lies between 7 and 15.5 pN. 

Are these forces consistent with atomic force microscopy (AFM) experiments for barnase (Best 
et ah, 2001)? First, let us caution that pulling a protein during mitochondrial import could be 
different than during AFM experiments which pull free proteins in solution. One mechanism for the 
difference is a drag force (Chauwin, Oster, and Glick, 1998). Another difference is that in AFM the 
displacement is imposed whereas during import the force is imposed. 

Best et al. (2001) conducted AFM pulling experiments on tandem repeats of barnase followed 
by the titin domain 127, and they extracted unfolding forces for the barnase units alone. Their 
data (circles) are plotted in Fig. 7 along side the data (error bars) of Table 1. The dependence 
of pulling rate k p on applied force F in the many AFM experiments to date follows Bell’s law 
(Evans and Ritchie, 1999): 

fc p = k° p e F ^ kT , (10) 

where kT is the product of Boltzmann’s constant and temperature, A:® and Az u are constants fit to 
the data. The quantity Aj: u has the interpretation of a potential-well-width for unfolding. Bell’s law 
plots as a straight line in a semi-log graph like figrefBest. In their experiments on the unfolding of 
titin domains, Carrion- Vasquez et al. (1999) observed that the zero force extrapolation (using Bell s 
law) of their AFM rates coincided with the zero denaturant extrapolation of chemical denaturation 
rates obtained for isolated domains of titin. This result suggests that extrapolation of the AFM data 
to smaller forces using Bell’s law might be valid. If the presently predicted forces are consistent with 
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Figure 7: Force vs. rate of pulling (AFM) or rate of import, o : AFM experiments (Best et al., 
2001); error bars: protein import data from Table 1. The five error bars from left to right correspond 
sequentially to the entries ir Table 1 


forces measured in the barnase AFM experiment, and both are consistent with the same Bell’s law, 
they should lie on the same ine in Fig. 7. This is difficult to conclude given the scatter in the data. 
All we can say at present is that the two are not inconsistent. Of course, this is not saying a whole 
lot given that a prediction of zero force would also appear to not be inconsistent when plotted on 
Fig. 7. 

A Net number of positive charges within the IM pore 

The net number of positive charges lying within the IM pore of the experiment (Huang, Ratliff, 
and Matouschek, 2002) is inferred as follows. First, from the pH of the experiment (= 7.4), we 
determined (from the Henderson-Hasselbach equation) that the ionization fraction is .998 or better 
for all the acidic and basic amino acids in the presequence. Since the pK of the N-terminus is 
uncertain (6.8 < pK < 8), its ionized fraction could range from .2 to .8 (with a positive charge); for 
simplicity we considered it as uncharged. For presequence (35), the experiment reported membrane 
potential catalyzed unfolding when positive charges were introduced at positions 15 or 16 but not 
when they were introduced at positions 18 or 19. Hence the last residue lying within the IM pore is 
either 16 or 17. Examination of the amino acid sequence then gives 3 net positive charges as lying 
within the IM pore. The number of charges for presequences (35; A16K) and (35; E15L) then follows 
naturally. Presequence (35+5) allows five more residues to occupy the IM pore and these add one 
more positive charge. The presequence (65) will have positions 46 or 47 (30 T either 16 or 17) as 
the last one lying in the IM pore. These consist of 6 or 7 net positive charges. However, taking 
the length of the IM pore to be between 65 and 90 A, we infer that between 20 and 28 amino acids 
(consisting of 4 to 6 net positive charges) will stick out at the matrix end of the IM pore. Hence we 
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have between 0 and 3 (6 to 7 minus 4 to 6) net positive charges. 
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